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We consider the evolution of the ground state in the Two Higgs Doublet Model during cooling 
down of the Universe after the Big Bang. Different regions in the space of free parameters of this 
model correspond to different sequences of thermal phase transitions. We discuss different paths 
of thermal evolution and corresponding evolution of physical properties of the system for different 
modern values of the parameters. 



The Standard Model relics on the Higgs mecha- 
nism of the electroweak symmetry breaking (EWSB). 
Its simplest realization (minimal SM) is based on a 
single weak isodoublet of scalar fields, which couples 
to the gauge and matter fields and self-interact via 
the quartic potential. The Two Higgs Doublet Model 
(2HDM) presents the simplest extension of the mini- 
mal scheme of EWSB, it contains two such doublets 
[l[ and a number of free parameters, for a review see 
e.g. 0- Many new phenomena can take place in 
this model in different regions of the space of these 
parameters. 

When describing the properties of 2HDM, one first 
finds the minimum of the Higgs potential and de- 
scribes its symmetry properties, then calculates the 
masses and the interaction of the gauge bosons and 
of the physical Higgs bosons, and then studies the 
Yukawa sector of the model. The phenomenological 
richness of the 2HDM is contained in particular in 
the fact that the Higgs potential can support minima 
(or, in general, extrema) of different nature. We will 
refer to them as the phases of 2HDM. Upon contin- 
uous change of free parameters of the potential, the 
exact position of the minimum shifts, but the sym- 
metry properties at the minimum remain unchanged, 
so that we still stay in a given phase. There are spe- 
cial points in the space of free parameters, at which 
this variation of parameters changes the nature of the 
global minimum, so that if the system crosses this 
point, a phase transition occurs. The phase tran- 
sition is reflected in a non-analytic behavior of the 
properties of the vacuum: vacuum expectation val- 
ues, masses, etc. 

At any given set of parameters, several extrema of 
the Higgs potential can coexist. Normally, we will 
be interested in the properties of the vacuum, i.e. 
the deepest minimum of the potential, however some- 
times it will be useful to trace the evolution of other 
extrema too. In particular, we will often calculate the 
difference in the potential depth between the vacuum 
(which we dub "the vacuum energy") and a higher 
lying extremum ("the extremum energy"). 

This richness can become even more manifest at 



non-zero temperature. The effective parameters of 
the model evolve with temperature, which can lead 
to thermal phase transitions, with important cosmo- 
logical implications. 

This issue has been studied in a number of papers, 
@) 0, H, S EH EH • All these papers focus on some 
more or less simple and specific variants of 2HDM 
and study them within thermal field theory at differ- 
ent levels of sophistication. In the present work we 
complement those papers by carrying out a system- 
atic analysis of all possibilities offered by 2HDM at 
finite temperature. We are interested in three basic 
questions (for an earlier discussion of these questions 
see [3): 

1. which possible thermal sequences of phases are 
allowed in 2HDM? 

2. what physical phenomena can accompany ther- 
mal evolution of the system? 

3. which modern, i.e. low-temperature, values of 
the parameters of the potential correspond to 
each phase sequence? 

Although we limit ourselves to the first non-trivial 
temperature corrections, we find novel phenomena in 
this approximation, which were overlooked in previ- 
ous works. Examples include the possibility of having 
an intermediate charge-breaking vacuum or new re- 
lations between the first-order phase transitions and 
the symmetries of the potential. 

In our analysis we want to stay as general as pos- 
sible and yet rely on explicit algebraic calculations. 
On the one hand, it is known that the most gen- 
eral 2HDM is intractable with straightforward alge- 
bra, so other methods must be used to gain some 
insight into its properties. On the other hand, it has 
been recently proved [l3[ that the entire list of possi- 
bilities for the thermal evolution of the most general 
2HDM is well preserved in a specific case of a soft Z<z- 
violating 2HDM, which we think is a more plausible 
candidate for the description of nature than the most 
general 2HDM. But in this case, analytical calcula- 
tions are possible, see e.g. [H, and we will use them 
to analyze all the sequences of phase transitions and 
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typical changes of physical properties of the system 
in detail. 

Let us also outline from the start the method that 
we will use. 

First, we limit ourselves to the first approximation 
of the thermal perturbation theory in the high tem- 
perature limit and to the tree approximation for the 
potential. This discussion is supposed to describe 
accurately thermal evolution far from the transition 
points and the phase sequences at sufficiently high 
temperature. At low temperatures, at the end of 
cooling down, the phase evolution of the Universe 
could have been more intricate than our discussion 
suggests. Still, we hope that our discussion is of some 
value even in this region as it provides a minimal list 
of phenomena possible in thermal 2HDM. 

The important advantage of this approximation is 
that the Gibbs potential has the same generic form as 
the zero-temperature Higgs potential, with only mass 
terms depending on temperature. This will allow us 
to describe thermal evolution of the Universe via a 
point moving in the phase diagram, whose structure 
we know from the zero-temperature analysis. 

Second, we will assume thermodynamical equilib- 
rium in each moment. This approximation is jus- 
tified when kinetic phenomena proceed fast enough 
and are not affected by the relatively slow expansion 
of the Universe. Near the second order phase tran- 
sition points, when the thermodynamical processes 
become slow, this approximation might break down. 

The structure of the paper is the following. In Sec- 
tion 1 we list some features of the 2HDM Lagrangian 
and present arguments in favor of the softly Z 2 - 
violating Lagrangian for description of reality. Sec- 
tion 2 is devoted to description of the temperature 
evolution of the parameters of the Lagrangian. In 
section 3 we briefly review all possible types of ex- 
trema of the potential and describe the phases of 
2HDM, following mainly 0, 0. In sections 4-6 
we consider, in the tree approximation, main phys- 
ical properties of all the possible phases in 2HDM, 
with a special attention to the case of explicitly CP- 
conserving Higgs potential, and we discuss possible 
ways the phase transitions could have occurred dur- 
ing cooling down of the Universe. Section 7 describes 
the emergent picture in general. 



I. LAGRANGIAN 

The electroweak symmetry breaking (EWSB) via 
the Higgs mechanism is described with the La- 
grangian 

C = C s g f +C H +C Y ; C h =T-V h . (1) 

Here, C s g f describes the SU{2) x U{1) Standard 
Model interaction of gauge bosons and fermions, 



whose form is independent of the realization of the 
Higgs sector; the Higgs scalar Lagrangian Ch con- 
tains the kinetic term T and the potential Vh > and Cy 
describes the Yukawa interactions of fermions with 
the Higgs scalars. 

In the minimal Standard Model (SM) one scalar 
isodoublet with hypercharge Y = 1 is imple- 
mented. Here the kinetic term has the standard 
form T = [D^cj))^ D^if), and the Higgs potential is 
V = 4>) 2 /2 — m 2 0*0/2. The minimum of V gives 
the vacuum expectation value v via ((f)) = u/yz = 
yJm 2 /2\. In this model there is one physical Higgs 
boson remaining; its coupling constants to the gauge 
bosons are expressed via their masses. The Yukawa 
interaction has form Cy = 9^ M Ql4"1r + h.c. with 
gf™ = V2m f /v. 

In the two-Higgs-doublct model the Higgs potential 
can have a rather complicated form. To describe it 
in a concise way, it is useful to introduce isoscalar 
bilinear combinations of field operators 



Xi = 4> 



x 3 = <j> 
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X2 
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— x 3 



(2) 



The Higgs potential of the most general 2HDM is 
conventionally parameterized as 

V = V -V 2 (x i ) + V i (x i ); 

V 2 {Xi) = MiXi EE 



- [mi 1 a;i+r7i22 ;z; 2 + (m 2 2 x 3 + h.c.)] , 
Vi(xi) = AijXiXj/2 ee 



(3) 



Xixl+X 2 xl 



+ X 3 Xix 2 + X i x 3 xl + 



+ X 6 xix 3 + \ 7 x 2 x 3 + h.c. 



At the first glance, this potential contains 14 free pa- 
rameters: real m 2 ^, m 2 ^, Ai, A2, A3, A4 and complex 
m\ 2 , A5, A6, A7. However the same physical content 
is described by a Lagrangian that can be obtained 
from ([3]) by a general rotation in the (0i, 02 ) space. 
This is known as the reparameterization symmetry. 
Therefore, physical observables depend not on all 14 
parameters but on a lower number of their combina- 
tions (cf. discussion in [l^| ) . An important particular 
case of reparameterization transformation is given by 
the rephasing transformation: 



{01 -> 01, 02 



9 2 C 
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A, 



A 5 e , 



etc. 



(4) 



Kinetic term. In the general case the renormaliz- 
ability requires that the kinetic term T is not diagonal 
but includes a mixed term II 611 



T = D, 



' 2 D^i 



(5) 
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This kinetic mixing can be eliminated by rotation and 
rcnormalization of scalar fields. However, in the gen- 
eral case parameters of this transformation vary with 
the renormalization scale even at small distances [l6| . 
This unattractive feature is absent in the case of ex- 
act or softly broken Z 2 symmetry, be it explicit or 
hidden. 

The Z 2 symmetry, exact and broken. At 

m i2 = 0, Ag = A7 = 0, >f = our system has a Z 2 
symmetry, i.e. it is invariant under transformations 

01 — ► <(>i , 2 — > ~4>2 Or 4>i — > -01 , 2 — > 02 ■ 

(6) 

In order for this property to survive through the per- 
turbation series, the Yukawa interactions must con- 
nect each right-handed fermion to only one scalar 
field 0i or 02 (Models I or II for Yukawa sector, see 
[H for details of the definitions). 

By definition, at A 6 = A7 = 0, h = 0, m 2 2 ^ 
the Z 2 symmetry is softly violated. In this case the 
parameters of the Lagrangian are determined unam- 
biguously up to the phase rotations of A5 and mf 2 
that keep Xtjn^ fixed. Models I or II for the Yukawa 
sector keep this property through the perturbations 
series. 

Hard violation of the Z 2 symmetry requires that 
at least one of quantities Xq or A7 or x differs from 
zero in each basis of the Higgs doublets. In this 
case for renormalizability all these coefficients appear 
via counter-terms, and the general Yukawa interac- 
tions couple each right-handed fermion to both scalar 
fields, 0i and 02 (known as Model III for the Yukawa 
sector). 

The general rotation of the Lagrangian of a softly 
i^-violating system in the (0i, 02) space induces 
nonzero Xq and A7, which makes the Lagrangian look 
as if hard violation of Z 2 were present. However, spe- 
cific relations arising among coefficients of the poten- 
tial prevent hard Z 2 violation (this situation is known 
as "hidden softly ^-violating case" ) . We do not use 
such a form below. 

Explicitly CP-conserving form of the poten- 
tial is the one that obeys the symmetry 0j «-> tf>\. It 
is provided by the condition that all coefficients in 
the Higgs Lagrangian ©, ([5]) are real. The general 
rotation in the (0i, 02) space makes this property 
hidden. 

Natural choice. It is often the case that a theory 
containing more than one field with identical quan- 
tum numbers arises in the low-energy limit of a more 
fundamental theory, in which these fields were com- 
ponents of a single multicomponent object. Such a 
theory often contains a higher symmetry, which might 
be not obvious at low energies. So, the fields, which 
had different values of the new quantum number as- 
sociated with this symmetry, should not mix at small 
distances, while large distances symmetry-breaking 
mixing terms can appear. If this argument is to be 



taken seriously for 2HDM, that is, if the quartic inter- 
action and the kinetic terms are supposed to reflect 
some new primary symmetry, then we should limit 
ourselves only with softly Z2- violating models, be- 
cause they allow mixing between 0i and 02 only at 
large distances. We think, therefore, that softly Z 2 - 
violating 2HDM is a natural, but sufficiently general 
choice to be considered [l6j . 

That is why we consider in detail evolution of the 
vacuum state for the Lagrangian with softly broken 
i^-symmetry in the explicitly CP conserving case 1 



V 2 (x. l ) = - 
Vi(xi) 



m n xi-r- 



m\ 2 X2- 



'12 



(353 + 4)) 



X 1 x 2 1 +X 2 xl f 
hA 3 a;ia;2-l-A4a;3a;^-|- 



A 5 (x| 



J2\ 



0. 



_ " (7a) 
According to [13j , this is a representative case of the 
most general 2HDM, in the sense that all phases and 
phase transitions that can happen in the most general 
case can be mapped to phases and transitions in this 
model. In particular, the results for the classification 
of regions in the parameter space and sequences of 
phase transitions derived below coincide with those 
obtained in the general case [l3| . 

Although we focus in this work on the scalar sector 
only, we note that softly broken ^-symmetric poten- 
tial is inherently stable in perturbative theory only if 
Yukawa sector has a form of Model I or II. Below we 
assume that is valid. 

It will prove useful to denote A2/A1 = fc 4 and in- 
troduce a special parametrization for the V2 term 



m\ x = m 2 (l — 5) , TO22 — k 2 m 2 (l + 5) , 



m\ 2 = \ikm 2 



k d ^ ^xT/xT 



(7b) 



so that we switch from the triple m\^, m 22 , m\ 2 to 
the triple m 2 , 8, /1. Note also that the type of the 
minimum realized at given m 2 j depends not on their 
absolute values, but on their ratios. Therefore, the 
parameter m 2 does not appear in the classification of 
phases, and we will need only to focus on the (fi, 5) 
plane below. 

To make some equations shorter, we also introduce 



1 The potential with hidden soft Z2 symmetry can be transformed 
to the potential of form Q by suitable a rotation. Besides, if 
CP-conserving extremum, with no scalar-pseudoscalar mixing, 
exists, then there exists a basis in the (</>i, $2) space, in which 
the potential has explicitly CP conserving form, with real A^, 
m ij' Hal' El- This is the reason why we consider the explicitly 
CP conserving potential l T7l . 



the following notation for certain combinations of A,;: 



A345 — A3 + A4 + A5, 

^345± = VA1A2 ± A345 



A345 — A3 + A4 — A5, 

-A-345± = VA1A2 i A345, 
A-3± = \AlA2 i ^3 • 



(8) 

The k symmetry. The quartic potential in ([7]), 
in addition to being .^-symmetric, is also symmetric 
under the transformation 

4>i^k4> 2 . (9) 

At 5 — this symmetry extends to the entire poten- 
tial. We call this property the k- symmetry. This k- 
symmetry is not the true symmetry of the Lagrangian 
since the kinetic term is not invariant under @, un- 
less k = 1. Loop corrections break this symmetry at 

Positivity constraints. To have a stable vacuum, 
the potential must be positive at large quasi-classical 
values of fields \4>k\ (positivity constraints) for an ar- 
bitrary direction in the (4>i, ^2) space. This translates 
into V4 > for all non-zero values of the fields, which 
places restrictions on possible values of A.;. For the 
potential© such restrictions have a simple form (see 
e.g. [13, El) 



Ai > 0, A 2 > 0, A 3 + > 0, 
A345+ > 0, A345+ > 0. 



(10) 



A geometric approach to the most general 
2HDM was developed in 14|, which helps gain some 



insight into the properties of the model without long 
cumbersome calculations. It is applicable to ther- 
mal evolution of the general 2HDM too, [13| , allowing 
one to see possible regions for existence of different 
types of vacuum and sequences of phase transitions. 
Since in this work we limit ourselves to the softly Z2- 
violating potential, we stick to a more conventional 
way of direct calculations, based on Q, to track down 
the properties of the system during thermal evolution. 



II. TEMPERATURE DEPENDENCE 

At finite temperature, the ground state of a sys- 
tem is given by the minimum of the Gibbs potential, 
which can be defined as 



V G = Tr 



(Ve~"/ T )/Tr(e-"/ T )^V + AV. (11) 



Corrections AV, to the first nontrivial approxima- 
tion, are given by diagrams of Fig. [Til Therefore, in 
this approximation the evolution of the system is de- 
scribed by the potential ([3]) with fixed quartic term V4 
and with the mass term V2 that evolves with temper- 
ature. This contribution is calculated with the Mat- 
subara diagram technique, giving a universal factor, 




FIG. 1: 1st correction to the Gibbs potential 

T 2 /12, for the loop integral at T 2 > m? [19]. The 
mass terms in the potential ([7]) can be then written 
as 



m 2 n (T) = m 2 n (0) - 2 Cl m 2 w , 
m\ 2 (T) = ml 2 (0) - 2k 2 c 2 m 2 w , 
m 2 2 (T)=mf 2 (0), 

T 2 



(12a) 



12m 2 



Here, the scalar loop contributions cf and the gauge 
boson loop contributions cf (Fig. |TTJ upper diagram) 
are 



3Ai + 2A 3 + A 4 



k 2 c\ 



s _ 3A 2 + 2A 3 + A 4 
c 2 - 2k 2 
(3<7 2 + 5' 2 ), 



(12b) 

with g and g' being the standard electroweak coupling 
constants. 

The fermion loop contributions c{ (Fig. |TTJ lower 
diagram) depend on the form of the Yukawa sector. 
For the Model II and Model I, the main contributions 
to these coefficients can be written in natural notation 

as 



c{(II) = yi ci(II) = -^g 2 b ; 
c f 1 (I) = l(g 2 +g 2 b ), c{(I) = 0. 



(12c) 



At given Ai these equations show that the curve of 
physical states is the straight ray in the 3-dimensional 
space (m 2 1; to 22 , w 2 2 ). This conclusion follows from 
the universal temperature dependence of the loop in- 
tegrals with 4>\ and <p2 circulating in loops in Fig. [TT| 

The corresponding equations for c\ in the most gen- 
eral case are presented in Appendix [A"l 



A. Ai space and (fj,, 5) plane 

The physical state of our system corresponds to 
some point in the Ai space, whose position is temper- 
ature independent, and a point on the (/x, 5) plane, 
which moves as temperature changes. As tempera- 
ture increases, this point traces a curve of physical 
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states on the (p, S) plane, starting from the "zero 
point", the point corresponding to the modern, zero- 
temperature state. Let us describe this curve explic- 
itly. 

Simple algebra allows us to express temperature- 
dependent parameters of the potential (|7b|) m(T), 
S(T), /i(T) via their "zero-point" (i.e. modern zero- 
temperature) values to, 5, /i: 

to 2 (T) = to 2 [1 - (c 2 + ciH , mCO = |U 



TO 2 (T) 



<5(T) = 



TO 2 (T) 



[<5 - (c 2 - ci)i 



These relations are transformed into equation 



<5 ( t) = aw {s 



P 



C2 ~ Ci 
C2 + Ci 



(13) 



(14) 



Eq. p4|) shows that the curve of physical states is a 
straight ray on the (/i, S) plane. These rays have also 
the following property: if extended back, beyond the 
"zero point", they all converge at a single point: 



Vo = (/i(T) = 0, S(T) = P) . 



(15) 



Of course, this point is physically unattainable, as it 
would correspond to T 2 — > — oo, however, its identi- 
fication is useful for visualization of possible thermal 
trajectories. Indeed, each ray of physical states can 
be constructed as starting form a certain "zero point" 
on the (6, /i) plane and moving strictly opposite to 
V . 

By virtue of the positivity conditions (fit))) the quan- 
tity C2 + C1 is positive 3 so that m 2 (T) decreases mono- 
tonically. Therefore, looking back into the past, with 
temperature increasing, we would see to 2 (T) decrease 
monotonically. As long as m 2 (T) stays positive, the 
coordinate n(T) along this ray grows from today's 
value \i > to infinity, with //(T)//x > 0. We will 
depict this evolution on the first sheet of the (fi, S) 
plane. At even higher temperatures, at to 2 (T) = 0, 
we pass onto the second sheet of the (/i, 6) plane, 
which corresponds to to 2 (T) < 0. The coordinate 
/j.(T)/fi in this ray is negative and grows with tem- 
perature rise from —00 to zero, while the point on 
this second sheet of the 6) plane moves strictly 
towards Vo- 

Below, in Figs. [5] and [SJ we will show these rays 
of physical states on the (p, S) plane in different sit- 
uations. We will see that sometimes these rays cross 



2 At k = 1 we have (neglecting fermions) ci = c\ => P = 0, i.e. 
5(T)/fi(T) = const. If additionally 6 = 0, we have S(T) = 0, 
i.e. the fc-symmetry is a symmetry of the Lagrangian, and 
temperature corrections do not break it. Note also that for a 

k 2 - 1 3VhM - 2A 3 - A 4 



purely scalar case P 



3 Indeed, c| 



k 2 + l Sv^iA^" + 2A 3 + A 4 ' 
oc 2A 3+ +A345+ + A345+ and cf > 0, c{ > 0. 



regions of different phases, which will have the form of 
ellipses or line segments. Directions of these rays, in- 
dicated by arrows, always correspond to temperature 
growth. Small dots, at which these rays start, indi- 
cate various possible "zero points". Note that once 
we have a ray, we have a freedom to place a "zero 
point" on it, so that any given ray can in fact corre- 
spond to several distinct zero-temperature situations. 
In all cases we select parameters of the potential in 
such a way that P < 0; rays for positive P can be 
obtained by appropriate mirror reflections. 



B. Sectors in \i space 

In accordance with the general geometric analy- 
sis of [3, the entire space of parameters Xi can be 
divided into sectors (in this approach each sector cor- 
responds to a relation among eigenvalues of some ma- 
trix), in which different types of vacua are allowed, 
but not necessarily realized: 



Sector I: A 345 _ > A 5 < A 4 + A 5 < 

Sector II: A 345 _ < A 3 _ < A 345 _ < 

Sector III: A 5 > A 4 A 5 > A 345 _ > 

Sector IV: A 5 < A 4 A 3 _ > A 4 + A 5 > 



(16) 



One can see that sectors do not overlap. Further in 
the text we will show how most of these constraints 
can be found from detailed analysis of extrema. 

In order to know which phase is realized at fixed 
set of Xi in each sector, one should look at specific 
regions on the (jj,, r))-plane. Thus, description of the 
phase diagram of 2HDM requires a two-step analysis 
— first select Xi, then study the (/x, 8) plane. 



C. Evolution of particle masses 

As the Universe was cooling down, the v.e.v.'s of 
the Higgs fields were changing. It may have resulted 
not only in variation of the absolute values of particle 
masses, but also in rearrangement of the particle mass 
spectrum, which can have interesting physical conse- 
quences. To this end, we discuss briefly evolution of 
masses of the gauge bosons and fermions, focusing on 
the Yukawa interaction in the form of Model II, Q. 

The mass of VF-boson is given by Mw = gv, 
while the quark masses are tud = 9dVi/\2, ra\j = 
guV2 1 v2) where g is electroweak interaction constant, 
go and gu are Yukawa constants and U = {u, c, t}, 
D = {d, s, b}. 

During thermal evolution all of v, V\, V2 change: 
v(T), vi(T), V2(T). We express high temperature 
values of particle masses via their modern values and 



G 



the ratio of v.e.v.'s to their modern values 



M W (T) = M w 



v(T) 



m D (l)=m D , mu(T) = mu- 



(17) 



Vn 



It follows that the ratios between masses of quarks 
of the same charge in different generations do not 
change, e.g. m t (T) : m c (T) : m u (T) = m t : m c : 
m u . However, relations between the quark masses of 
different charge can vary during evolution, 

m D (T) _ m D cos/3(T) mu(T) _ mu sin /3(T) 
M W (T) ~ Mw cos/3 ' M W (T) ~ ~Mw~ sin/3 ' 

mu(T) mu tan/3(T) 

mjj(T) mu, tan/3 

(18) 

We will show several graphs below that illustrate be- 
havior of the particle spectrum in various situations. 



III. EXTREMA OF POTENTIAL. MAIN 
TYPES 



The Higgs potential can have several extrema. The 
extremum with the lowest value of the energy, the 
global minimum of potential, realizes the vacuum 
state. The other extrema can be either saddle points 
or maxima or local minima of the potential. 

The extrema of the potential define the values 
2) of the fields 4>\,2 via equations: 

dv/dfa\^ =i<Pz) =0, dv/d4\^ { ^ =0. (19) 

These equations have the two main types of solutions: 

• electroweak symmetric (EWs) solution 

(4>i) = 0; 

• the electroweak symmetry violating 
(EWv) solutions with at least one (fa) =^ 0. 

In the next several Sections we will describe these 
phases and possible thermal phase transitions in the 
following sequence: 

• EWs solution and EWSB phase transitions; 

• classification of EWv solutions; 

• evolution of phases in the case when only 
EWSB phase transition takes place; 

• evolution of phases in the case of a first-order 
phase transition between two CP-conserving 
phases; 

• evolution of phases across CP- 
violating/restoring phase transition; 

• evolution of phases passing through the charge- 
breaking phase of 2HDM. 



IV. EWS PHASE. EWSB PHASE 
TRANSITION 

The EWs point (fa) — (fa) — is an extremum of 
the potential. In its vicinity only the mass term V2 
can be used for the analysis of the stability. 

The EWs point is a minimum of the potential if 

m\ x < > m 22 < an d > l TO 12| 2 ) 

(20a) 

i. e. within the circle lying on the second sheet of the 
(/1, S) plane 

C:/z 2 + r5 2 = l with m 2 (T) < => v(T)/fi < 0. 

(20b) 

The positivity constraints guarantee that this mini- 
mum realizes the vacuum 4 . 

Fig. [5] shows possible patterns of high-temperature 
evolution of physical states, which takes place on the 
second sheet of the (//, 6) plane. The rays on this 
second sheet are shown with solid lines; they are the 
high-temperature continuations of the rays from the 
first sheet shown as dotted lines (which are to be 
described later in Figs. [5] and |5J). As temperature 
rises, the point moves towards the special point Vq, 
reaching it at infinite temperature. If it intersects 
the circle C (/.t 2 + S 2 = 1) at some temperature, a 
standard EWSB phase transition takes place. 

As can be seen from (TT21 , the evolution of m\ x and 
7Ti22 i s m general different, and they change their signs 
at distinct temperatures. Thus, in contrast with the 
minimal SM, we can have an intermediate situation 
when one of m\ is positive, while the other is nega- 
tive. In this case EWs extremum is a saddle point. 

Fig. [3 shows several scenarios of how the EWSB 
phase transition could have taken place at different 
parameters of the potential. 

• If ci > and c 2 > 0, then \P\ < 1 and the 
point Vo (IT51) lies inside C. As temperature 
increases, the point will always cross the cir- 
cle exactly once on its way towards Vo- This 
happens no matter what the zero-temperature 
parameters are. This is shown in the left panel 

of Fig. m 

• If C1C2 < 0, Vo lies outside C. In this case the 
system will reside in the EW-violating phase 
even at extremely high temperature. With an 
appropriate choice of modern values of /z, 5, the 
ray of physical states can cross the circle C, but 
it does so twice, so that the EWs phase becomes 
an intermediate phase. This is illustrated in the 



4 At [i{T) I \i > 0, i.e. on the first sheet of (/t, S) plane, and inside 
the similar circle l !20bl , the EWs point is a local maximum 
of the potential. For all points outside this circle the EWs 
extremum is a saddle point. 
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FIG. 2: The second sheet of the (fi, 5) plane. Dotted rays are curves of physical states at the first sheet of this plane, which 
begin at the modern zero-temperature point (jig, <5o) or (fig, S' ). Solid lines are their high-temperature continuations on 
the second sheet. The left and right panels correspond to \P\ < 1 and \P\ > 1, respectively. At points where the solid 
lines intersect the circle C the system experiences EWSB phase transitions. 



right panel of Fig. [2j ray 2. Another possibility 
is that there is no phase transition at high tem- 
perature at all, see ray 1 in the same picture. 
In the following sections we discuss only the case with 
a single EWSB transition, \P\ < 1. The analysis of 
the other possible cases for m 2 (T) < can be con- 
ducted in a straightforward way, but it is decoupled 
from the discussion of the first sheet of the (/i, S) 
plane. 



ELECTROWEAK VIOLATING (EWV) 
EXTREMA 



For the EWv solutions we consider also the values 
of operators Xi at a specific extremum point, which we 
denote generically as N (we will omit such subscripts 
whenever it does not lead to confusion). 



2/1, n 
2/2, N 
U3,N 
2/3, TV 



(xi)n 

(X3)N 
(4)jV 



hy N (<i>2)n , 



Pl/N , 



V2) N 



1 (4>l)N 



At each extremum the quantities 2/1,2 are positive, 
and the Cauchy inequality holds for an important 
quantity Z: 



yi > , ?/2 > . 



V1V2 -y* 3 y 3 >0. (21) 



Thanks to the properties of homogeneous 
functions, we have at each extremum 
V 2 {(<f>i) N ) = -2V4({<t>i) N ). Therefore, the ex- 
tremum energy is 

£ff = V(((f>i) N ) = = V 2 {{4> l ) N )/2. 

(22) 

For each EWv extremum one can choose the z axis 
in the weak isospin space so that (<fii) = 

with real positive V\ (choose "neutral direction" ) . In 
this basis ((^2) has in general an arbitrary form. Then, 



after this choice the most general electroweak sym- 
metry violating solution of (|19p can be written in a 
form with real v\ and complex V2' 



with v\ 



= -( ° 

V2 \vi 

= \ v i\ = wcos/3, V2 



_ J_( u 

|x?2 |e z ^ = v sin /3e^ . 

(23) 

Without loss of generality one can consider only real 
positive u. 

For the most general potential ([3]) one can classify 
the EWv extrema according to their values of Z [i| : 

• If u 7^ 0, then Z > — charge-breaking 
extremum. 

• If u = 0, then Z = — neutral extremum 

with v 2 = 2(yi + y 2 ). 

A Higgs potential with arbitrary complex As^,?, 
m 2 2 (i.e. an explicitly CP violating potential) gen- 
erally possesses extrema, at which neutral scalars 
have no definite CP parity (in these states generally 
£ = arg«2 7^ 0). Such potentials in general do not 
offer any helpful classification of the extrema. 

In the case of explicitly CP conserving potentials, 
which we focus on in this work, one can further clas- 
sify neutral extrema as 

• CP conserving — CPc, 

• spontaneously CP violating — sCPv 

Since there can be up to two CP-conserving local 
minima with different properties, we will refer to var- 
ious phases as CPcl, CPc2, sCPv, charge-breaking. 

Two properties of 2HDM, proven in ref. (TiJ, are of 
much importance for future discussion: 



No more than two distinct local 
minima of potential can coexist. 



The minima that preserve some symmetry 
and that violate it cannot coexist. 



(24a) 



(24b) 



Property (|24ap means in particular that if two min- 
ima are degenerate, they realize the vacuum state, 
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FIG. 3: Parameters of the vacuum state for the case with a single phase transition EWs^CPc. Shown are the vacuum 
and extremum energies (upper left), tan/? (upper right), as well as v\, V2 and v in the lower row. 



and the other extrema arc not minima. 



It follows from the property ()24b|) that transitions 
between the states of different symmetry is realized 
via merging of these extrema. Therefore all parame- 
ters change in this transition continuously, that is the 
corresponding phase transition is of the second or- 
der (certainly, taking into account fluctuations might 
modify this conclusion). 



VI. THE EWV VACUA AND PHASE 
TRANSITIONS 



The equations for v.e.v.'s and extremum energies 
were obtained in Q for each type of extremum. They 
allow to obtain equations for the sectors in the param- 
eter space (|16l) . where certain types of vacuum are 
realized. Note that transitions among various EWv 
phases can take place only at m 2 (T) > 0, i.e. on 
the first sheet of the (/i, 6) plane, which are shown in 
Figures below. 

CPc extrema, general. The CPc extrema are re- 
alized in the entire space of parameters of the poten- 
tial. When searching for such an extremum, one can 
transform two cubic equations representing condi- 
tions (fl"9]) into relations for quantities v 2 — 2(y\ +1/2) 
and r = kt = fctan/3 = kyyzjyl. For potential 



these equations can be written as 



v 2 = m 2 (k 2 +T 2 ) ' '' + / ' 7 



A345T 2 + \A1A2 ' 
VA7A^ 4 + (A 3 45--£A 34 5 + )t 3 - (25) 



-(A 345 _ + <5A 345+ )t - = . 

An algebraic equation of fourth degree can have 0, 2 
or 4 real solutions (including accidental degeneracy). 
By construction, one should consider only real solu- 
tions of equation (f2"S"l) satisfying v 2 > 0. Omitting 
solutions with possible negative values of v 2 , one can 
state carefully that there could be up to 4 CPc ex- 
trema. 

The extremum energy is defined via solutions of 
these equations as 



&CPc — 



l k 2 (1-<5+/xt)[1 - 5 + 2^r + t 2 (1 + 5)} 



A 3 45T 2 + VA1A2 



(26) 

Let us now specify different sectors in the A^ space, 
and discuss the characteristic properties of thermal 
evolution of the system in each sector. 



A. Sector I. No phase transition except EWSB 

We start with the case, when there is no other 
phase transitions apart from the EWSB one. Such 
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FIG. 4: The same as in Fig. [3| but for another set of \i . 



a situation is possible for any A.; , within an appropri- 
ate region on the (/x, S) plane. Besides, in sector I 
in the Xi space defined by 

A 5 < , A 4 + A 5 < , A 345 - > , (27) 

this situation takes place for any "zero point" on the 
(fi, S) plane. These inequalities will be proved below, 
when we will show that in sectors II, III, and IV sev- 
eral phase transitions can happen, and then (|27|) will 
appear as a complement to those sectors. 

The only phase transition possible in this sector 
is EWSB transition between the EWs and CPc ex- 
trema. Figures [31 HI represent, for two different sets 
of Xi , the thermal evolution of phases and their prop- 
erties in these cases. 

Each set of graphs contains evolution of sev- 
eral physical quantities as functions of temperature- 
dependent parameter w oc T 2 . The upper left plot 
shows the evolution of the vacuum energy (thick 
curve) as well as the energies of other extrema (thin 
curves). The dashed line indicates the position of 
the EWSB transition. The upper right plot gives the 
evolution of tan/3. The lower row contains dependen- 
cies of v\ and i>2, which within Model II for Yukawa 
sector, give also the temperature dependence of the 
up- type quark and down- type quark masses. The last 
plot shows evolution of v, which is proportional to the 
mass of the gauge bosons. 

It is remarkable that even without phase transi- 
tions, thermal evolution experiences stages of very 
fast change of the quantities shown. Another re- 



markable observation is that sometimes masses of the 
quarks and gauge bosons can depend on temperature 
in a non-monotonic way. A possibility of rearrange- 
ment of the quark mass spectrum is clearly seen from 
the plots of vi and v-i. 

B. Sector II. First order phase transition 
between two CPc vacua 

A first order phase transition between two CP- 
conserving phases can appear only in the sector II 
in the Xi space, which is bounded by the following 
inequalities: 

A 345 _ < , A 3 _ < , A 345 - < . (28) 

We explain the origin of these inequalities in this sec- 
tion below. 

A phase transition between two phases of identi- 
cal symmetry implies that at some temperature these 
phases are degenerate. It was proved in [l4| that such 
a degeneracy always appears in the tree level 2HDM 
potential as a result of spontaneous breaking of an 
extra symmetry of the potential (and not necessarily 
of the whole Lagrangian) 5 . In our case this is the 
fc-symmetry, which we introduced before. One can 



5 The explicit form of the fc-symmetry changes under 
reparametrization into more complex linear transformation. 
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prove that there are no other equally generic symme- 
tries possible. 

• Degenerate case with fc-symmetry ([9|) is re- 
alized at 5 = in ([7]). In this case the equation for r 
([25]) breaks in two independent equations 



(r 2 - 1) (v/AiA^t 2 + 1) + A 345 -t) 



= 



which can be easily solved. 

A. The first of these equations yields two non- 
degenerate solutions A±: 



(A) : t 2 = 1 => t = t a = ±1 



^CPcA — 



m 4 k 2 (1±^) 2 



(29) 



A 



345+ 



B. The second equation yields two degenerate 
solutions B±: 



(B) : r 2 



A345- 



:T + 1 = => 



r = tb 4 



2^X1X2 



£cPcB ± 



W ( 1 



4 V2VaTa^ a 



345- 



345- 

(30) 
(31) 



These solutions show that fc-symmetric vacuum can 
exist only along the segment 



2y A1A2 



(32) 



in the (fi, <5)-plane (Fig. [5]). 

Let us calculate now the difference between ex- 
tremum energies (f3"Tj) and ([2"9"]l . Remembering that 
2VA1A2 = A345_ + A345+, one can obtain easily 



£cpcB ± —£-CPcA 



2 k 2 [A 3 45-(l±/l)±A 3 45+Mr 



4 2 A 345 _ A 345+ V Ai A 2 



From this it follows that the state B can have a lower 
energy than the state A, only if A 345 _ < 0, which 



gives the first condition in (j2"8|) together with useful 
condition for vacuum state > fits- > 0. 

Now, a direct calculation of masses of the charged 
and pseudoscalar Higgs bosons for solution B gives 

M 2 H± = -v 2 A345 - /2 , M\ = ~v 2 A 3 - 12 . 

These quantities must be positive at the minimum of 
the potential, which results in the second and third 
conditions of ([28]). The conditions (J28J) and (J32J) form 
necessary and sufficient conditions for realization of 
vacuum B±. 

For solutions B, it is useful to calculate separately 
v.e.v.'s for each field <bi\ 



21,2 



1 „ m 2 k 

v l = 



l + t 2 



u 2± 



l + t 2 VAiX 



O 2 ( 33 ) 
-.(1±R). 



Weak violation of the fc-symmetry (quasi- 
degenerated case). Phase transition. 

As temperature changes, the quantity S goes through 
zero, see ([T4"]) . It is useful to consider the case of small 
5. In this case corrections to yi and t are easily cal- 
culable from ([2"5"]) considered as solution ([3TJ]) with a 
small correction. We do not present here these equa- 
tions but show just the corresponding value of ex- 
ternum energy in the following way: 



m k 



41.2 



CPcB± 



[terms identical for B±— 



2^ 2 VXlAI- A 



345- 



A345- (M 2 VA1A2" + A345) 



■6 



(34) 

This equation shows that as the curve of physical 
states passes through (5 = 0, the phase with the lowest 
extremum energy switches from B + to £?_ (or vice 
versa) . That is a phase transition from the phase B± 
to the phase B^ with a jump Vi+ — > Vi—, V2+ — > U2— • 
Since the symmetries of these two phases are identical 
even at 8 ^ 0, theorems ([24]) allow for existence of two 
distinct minima of the potential at small enough 5: 
the vacuum and a meta-stable state. 
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FIG. 6: Evolution of the parameters of the vacuum state for sector II, ray II. Upper row: Extremum energy, tan (3, total 
vacuum expectation value v. Lower row: Vi and V2 shown separately. 



The discussed phase transition is a first order 
phase transition. The latent heat of this phase 
transition is given by the well known thermodynam- 
ical equation (see Cj in sect. [TTJ) 



= T- 



M4 



- T- 



M- 



dT 

- A 3 4 5 _ 



dT 



/i 2 VAiA 2 + A 



:S45 



5^0 

R 



VA1A2 



(ci - c 2 )w . 
(35) 

The case with non small 5 can be studied by 
directly solving eq. (f25]) . These solutions are direct 
continuation of above quasi-degenerated case. 

• The (/i, 8) plane. Fig. [5] represents the first 
sheet of the (fi, 6) plane for some specific sets 
of parameters from sector II. Several possible 
"zero point" states (i.e. different choices of zero- 
temperature /i and 8) are shown. The first order 
phase transition segment, described by eq. (f3"2")l , is 
marked by thick line. Possible evolution paths of 
physical states are presented by rays, arrows indicat- 
ing the direction of temperature growth. 

Small dots on these rays labeled 1,2, and 3 cor- 
respond to possible "zero points". The shaded area 
covers all "zero points" that would lead to crossing of 
the thick segment at some temperature, which would 
force the system to go through a first order phase 
transition (e.g. point 2). In this case the phase evo- 
lution during cooling down of the Universe was 
EWs -> CPcl -> CPc2, 



that is, the Universe experienced two phase tran- 
sitions: EWSB, of the second order, and CPcl — > 
CPc2, of the first order). 

If today's state corresponds to the "zero point" 3 
lying on ray II, or if it is located anywhere on ray I, 
then the Universe has experienced only one standard 
EWSB phase transition during its cooling down 6 : 

EWs -► CPc. 
This case was described in the previous subsection. 

• Evolution of the parameters of vacuum 
state. A typical evolution of the parameters of the 
vacuum state for ray I has already been shown above 
in Figs. Figs. |H Figs. [H [7] represent evolution of 
the parameters of the vacuum for ray II starting from 
the "zero point" 2. The same graphs can be also used 
to track the evolution of the same ray, but starting 
from another "zero point" 3. In all these Figures 
we show evolution of the parameters of the vacuum 
state as functions of the parameter w (p~2|) , the tem- 
perature increasing from left to right. The evolution 
of the physical parameters of the true vacuum state 
is shown everywhere by thick lines, while the thin 
lines correspond to the other extrema and are shown 
just for comparison. The phase transition points are 



6 Note that at k = 1 we have c\ = C2 => Vo = (0, 0). The 
ray of physical states never crosses the line segment of the first 
order phase transitions <5 = 0, so, we have no phase transition. 
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FIG. 7: The same as in Fig. [5| but for another set of Xi. 



indicated by small stars. 

In all pictures the upper left plots represent ther- 
mal evolution of the extrema energies shown with 
thin lines. The vacuum state is drawn with thick 
line. The upper middle plots show the evolution of 
tan/3 = V2/V1, while the upper right plots gives the 
temperature dependence of the main v.e.v. value 
v = \Jv\ + v\ = Mw/g- The lower row shows the 
two v.e.v. 's Vi and V2 separately. Within Model II 
for the Yukawa sector v\ oc m u (rrit), i>2 oc md {nib)- 

The first order nature of the phase transition is ev- 
ident on these plots, as the quantities v, Vi and tan/3 
jump at the transition point. These changes can be 
large in their magnitude. One can see that for the 
different sets of Xi, distinct variants of thermal evo- 
lution can be realized. For the cases shown in Fig. [5] 
(Fig. [7]) the following changes take place at the phase 
transition, respectively: 

• tan/3 jumps up (down), 

• main v.e.v. v as well as the masses of the gauge 
bosons My/, Mz jump down (up), so that the to- 
day's value of v can be either larger or smaller than 
before the phase transition. Similar behavior can be 
demonstrated by Vi separately. In particular, Fig. [JJ 
displays a peculiar My/ oc v behavior as tempera- 
ture increases: at first it monotonically decreases, 
then jumps up at the first order phase transition, and 
then goes down again and turns zero at the point of 
EWSB. 



C. Sector III. Spontaneously CP violating 
(sCPv) minimum 

The spontaneously CP violating vacuum can exist 
only in sector III, which is defined by the inequali- 
ties 

A 5 > 0, A 5 >A 4 , A 345 ->0. (36) 

We again will explain the origin of these limitations 
in this section. 

Since the potential is CP symmetric, the sCPv ex- 
ternum is doubly degenerated in the sign of the CP 
violating phase £. Because of the property ([24]) . if 
a sCPv minimum exists, it is necessarily the global 
one, i.e. the vacuum (for more details see [H). 

To describe a sCPv extremum, it is useful to ex- 
press the potential (UJ via quantities y\, yi and 
£. This gives a second-degree polynomial in cos£. 
Therefore the extremum condition can be solved first 
for cos£. It has two solutions, the CPc solution 
sin£ = and sCPv solution with £ 7^ 0. We dis- 
cuss here only the sCPv solution. The coefficient in 
front of the cos 2 £ term in the potential is propor- 
tional to A5, so if we want sCPv to be the minimum, 
we must require than A5 > 0, which yields the first in- 
equality of (jnSJ) ■ Besides, a direct calculation of the 
charged Higgs boson mass for the sCPv phase, [1], 
gives Mfj± = (A5 — A4)w 2 /2, whose positivity leads 
to the second inequality in (|36|) . 

Inserting cos£ obtained back into the extremum 
condition yields a second-degree polynomial in yi.2- 
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FIG. 8: The (n, 5) plane for sector III. 




FIG. 9: Evolution of the parameters of the vacuum state though the double CP-violating and CP-restoring phase 
transition. Upper row: Vacuum energy evolution (left), the insert zooms at the region between the two phase transitions; 
tan/3 (center); v (right). Lower row: vi (left); V2 (center); sin 2 £ (right). 



The minimum conditions for yi take form of two si- 
multaneous linear equations, which have unambigu- 
ous solutions for any values of the parameters (with- 
out assumption on Z2 symmetry). Finally 



1)2- 



k 2 m 2 



m 



1 



A345+ 
1 



A345- 
s 



cos £ = Co 



A345+ A_345_ 

fikm 2 



m 4 k 2 



sCPv 



4A 5v / yiy2 
1 S 2 



_A-345+ A 3 4 5 _ 

Equation (J37cj) has two solutions 

= arccosco , £- = — 



JJ_ 

2A 5 



(37a) 

(37b) 
(37c) 

(37d) 
(38) 



The vacuum energy E s c p v (|37d|) is the same for both 
these solutions: the vacuum state is doubly degener- 
ate in the sign of the CP violated phase £. 

Simple algebra allows one to rewrite the condition 
cos 2 £ < 1 as the condition that the sCPv state lies 
only inside a specific ellipse in the (//, S) plane: 



[i 2 S 2 

— — I = i 

b\ + b 2 



h = 



2A fl 



A. 



62 = 



A 



345- 



345+ 



A 



(39) 



345+ 



A geometrical analysis shows that the necessary con- 
ditions for the sCPv vacuum contains the condition 
^345- > 0, which is the third inequality in (f31)|) . [131 ]. 

• The (fi, S) plane. Fig. [H] shows the first sheet 
of the (//, 5) plane for different sets of Xi from sector 
III. The states of the sCPv vacuum are located within 
the gray ellipse, described by eq. (|39|) . The left plot 
represents a case when the point Vq lies outside the 
ellipse, while the right plot refers to a situation when 
the point Vq is inside the ellipse. 

Possible evolution of physical states is again 
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FIG. 10: The same as in Fig.^for another set of parameters. 



presented here by rays, with arrows indicating the 
direction of the temperature growth. The hatched 
area covers all the modern values of (/i, 5), such 
as point 2 on ray II in left plot, for which thermal 
evolution in the past went through the ellipse of 
sCPv vacua with two phase transitions. For these 
states phase evolution during cooling down of the 
Universe was 

EWs — ► CPc -> sCPv -> CPc 
(three phase transitions of the second order). 

Today's states such as 3 describe a sCPv vacuum 
(for both plots in Fig.[8|). For these states the phase 
evolution during cooling down was 

EWs -► CPc -► sCPv 
(two phase transitions of the second order). 

If today's state corresponds to the "zero point" 
4 on ray II, or if it is located anywhere on ray I, 
then the Universe during its thermal evolution had 
just one standard EWSB phase transition, the case 
already described above. 

• Evolution of parameters of the vacuum 

state. Thermal evolution of the parameters of the 
vacuum for ray I is the same as in Figs. [3l [4] We 
discuss here the evolution of the parameters of the 
vacuum state along ray II of Fig. [HJ which is shown in 
Figs. 191 HHl and which reflects a series of second order 



phase transitions CPcl^sCPv^CPc2 (the phases 
CPcl and CPc2 can be either different or identical). 
Notation and description of these plots are the same 
as for the previous case Fig. [6] An essentially new 
picture, shown last in the lower row, represents the 
behavior of the order parameter for sCPv case, sin 2 £, 
which is, of course, non-zero only for the sCPv phase. 

It is interesting to note that for the set of chosen 
in Fig. [9j the energy difference between the vacuum 
and the second deepest extremum is small and prac- 
tically invisible on the main graph. It makes fluctua- 
tions from the vacuum sCPv state into the low-lying 
CPc states (saddle points) not that much suppressed. 

The order parameter sin 2 £ exhibits a similar be- 
havior in both cases. Thermal evolution of the other 
parameters makes it evident that we deal with second 
order phase transitions. The curves are continuous 
but their slopes experience jumps at the transition 
points (indicated by small stars), similar to e.g. sus- 
ceptibilities in condensed matter physics. 

All other features of evolution of the physical pa- 
rameters can be different as in previous cases. V.e.v.'s 
can either decrease monotonically or increase at some 
temperature. In particular, for the case of Fig. [9] 
with going to the past in the CPc phase the mass 
Mw oc v decreases (together with mass scale parame- 
ter m 2 (T)) but in the sCPv phase it grows strong with 
increase of temperature. After second phase transi- 
tion all masses become decrease monotonically to the 
past, starting from new level, heavier than modern 
(for Mw and m u ). 
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FIG. 11: The (/x, 6) plane for sector IV 



D. Sector IV. The charge-breaking vacuum 
(with u ^ 0) 

The charge-breaking vacuum can exist only in sec- 
tor IV, which is defined by the following inequalities 
(fore more details see Q): 

A 4 ± A 5 > , As- > . _ (40) 
In a charged vacuum, it is not possible to split 
the Higgs boson and the gauge boson mass matri- 
ces into the neutral and charged sectors, the inter- 
action of gauge bosons with fermions will not pre- 
serve electric charge, photon becomes massive, etc. 
[2ll ]. This is the reason why such vacuum is called 
charge-breaking. Certainly, this phase cannot de- 
scribe today's Universe, however the vacuum could 
have evolved through this phase in the past. 

In the problem of finding the charged extremum 
one can treat variables yi as independent, so that 
the position of this extremum is given by the system 
of linear equations, written in variables yt, dV/dyi = 
=> Ay j/j c/j = Mi. If Ay is not singular, a solution 
to this system always exists and is unique (but it 
might violate conditions (|2ip). 

Decomposition (|B1I) shows that if such an ex- 
tremum realizes the minimum of the potential, this 
minimum is the global one (vacuum) . This decompo- 
sition shows that the charged extremum can be min- 
imum if only V&{xi) > for all classical Xi regardless 
of conditions (I21|) (the necessary condition). In par- 
ticular, for the potential (J7]) it means that the con- 
ditions (fTU|) are to be accompanied with inequalities 
(|4"0|) . For this potential the extremization problem is 
solved easily, and after simple algebra the solutions 
can be written [4| in form 



Vi = 

V2 



m 
~2~ 



cext 
c ch - 



2/3 

m 4 k 2 



1 

A 3 + 
1 

A 3+ A 3 
m 2 fc/i 
= 2(A 4 + A 5 ) ; 

1 : JP_ 
A3 j 



A 3 - 
S 



(41) 



A 3 _ A 4 + A 5 



which is necessary for the existence of a charged min- 
imum (vacuum), defines the interior of an ellipse in 
the (/i, 5) plane: 



62 1 
— < 1, 



where ai 



Simple algebra shows that condition Z > (|2"Tj) , 



A 4 + A 5 A 3 _ 
a 2 A 3+ A 3+ 

(42) 

At m 2 (T) > and within this ellipse we also have 
yi > 0, Tj2 > 0, therefore the charged minimum exists 
and realizes the vacuum. If /x, S point lies outside this 
ellipse, the minimum is neutral. 

• The (/i, S) plane. Fig. QTj represents the first 
sheet of the (fi, S) plane for different sets of A,; from 
sector IV. The states of charged vacuum are located 
within the grey ellipse described by eq. f|42[) . The left 
and right plots represent the case when the point Vo 
is outside or inside the ellipse, respectively. One can 
see that the situation is completely the same as in 
the sCPv case (Fig. |8j , with only difference that the 
ellipse here represents the charged vacuum instead of 
sCPv. 

The shaded area covers all "zero points" on the 
(/1, S) plane, which crossed in the past the ellipse 
of the charge-breaking vacua (points 2 at rays II 
in the left plot only). Such trajectories describe 
the following phase evolution the during Universe 
cooling down 

EWs -> CPc -> charged -> CPc 
(three phase transitions of the second order). 

We obviously disregard the cases when the "zero 
point" lies inside the charged vacuum ellipse. 

In the case of "zero point" 4 on ray II, if the "zero 
point" is located on ray I, during cooling down, the 
Universe undergoes one standard EWSB phase tran- 
sition EWs -> CPc (see sect. IvTa"|) . 

• Evolution of the parameters of the vacuum 
state. A typical evolution of the vacuum parameters 
for ray I was already shown in Figs. 02 Figs.Hl Fig. [12] 
represents evolution of the parameters of the vacuum 
corresponding to ray II with "zero point" 2; evolution 
from "zero point" 3 can also be read off these plots. 
The notation and description of figures here are the 
same as for previous cases, e.g. in Fig. [S] The differ- 
ence is that we do not show evolution of v.e.v.'s in the 
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FIG. 12: Evolution of the £ parameter and other physical quantities during transition through the charge-breaking 



charged vacuum phase, since these quantities have no 
clear physical sense in this phase. Instead of that, in 
the right bottom plot we present temperature depen- 
dence for the order parameter in the charged vacuum 
phase, which we choose as £ = y/ Zj (1/32/3) = u / v 2, 
i.e. the relative value of charge symmetry violating 
component of v.e.v. (fa)- When the ray of physical 
states crosses the boundary of the ellipse P2")l . this 
parameter starts growing from zero (Fig. [T2")), which 
is a typical behavior of an order parameter of a second 
order phase transition. 

The evolution of the v.e.v. 's near the charged vac- 
uum phase demonstrates the same features as for the 
series of phase transitions through the sCPv vacuum 
phase (Figs. HO}. 



E. What if potential is not explicitly CP 
conserving? 

In Section U we presented arguments in favor of the 
softly Z2 violating form of the Higgs potential. In 
addition, to make discussion more specific, we lim- 
ited ourselves to explicitly CP-conserving potentials 
(0 . Here we briefly discuss what changes if the latter 
condition is lifted. 

In this case a softly Zi violating potential has form 
([7]) with complex A5 and m\ 2 - The rephasing trans- 
formation (j4|) allows one to make A5 real and consider 
only ji complex. This approach allows one to use the 



same subdivisions of space into sectors as we used 
in this Section. 

According to [l5|, [2f3| . in this case the CPc phase 
does not exist. All phase states violate the CP sym- 
metry, and we will denote them as CPvl, CPv2, 
etc. Straight rays on the (/i, 8) plane are replaced 
by straight rays in the 3-dimensional (i?e/i, 7m/z, 8) 
space. 

The situation in sector IV is analyzed in the same 
way as before. With equations from [4], the condi- 
tion for existence of the charged vacuum within the 
ellipse (|42l) is transformed into a similar equation for 
3-ellipsoid. Therefore, the phase evolution is the same 
as discussed in sect. IVIDI 

For other sectors a detailed study with explicit gen- 
eral equations for physical quantities and their behav- 
ior in different phases becomes very complicated. In 
these cases the geometrical analysis of [13[ allows one 
to obtain general picture of phase transformations. 
In sector II and in the case of no other phase transi- 
tions except EWSB, we obtain roughly the same pic- 
ture as above, sect . I VI XI I VI B I with a natural change 
of the labels: CPc^CPv, CPvl, CPv2. In sector III 
the picture changes, as the doubly degenerated sCPv 
states disappear. The ray of physical states goes now 
in 3D and is inclined towards the plane, where an 
ellipse of doubly degenerate states similar to sCPv 
exists. The ray can traverse this plane, resulting in 
a first order phase transition, like the one in our sec- 
tor III, with a nonzero specific heat that depends on 
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the slope of this ray. Nevertheless, the region of all 
modern values of (fx, S), which would go through the 
ellipse at high temperature, can be described in an 
essentially similar way as the shaded area in Fig. [5J 

VII. DISCUSSION 

We considered here the phases and thermal phase 
transitions in the 2HDM. We limited ourselves with 
the case of softly broken Z2 symmetry, presenting ar- 
guments in favor of its realization in Nature. Besides, 
in what concerns sequences of phase transitions this 
case is representative of the most general situation 
[l3j . obtained in tree approximation. We obtained a 
rich picture of possible phase states and phase tran- 
sitions. Taking into account finite mass and higher 
order effects will modify the simple temperature de- 
pendence of effective potential parameters discussed 
in sect. QlJ We believe that with these modifications 
the picture of possible phase states and phase transi- 
tions can only become richer. In particular, some sec- 
ond order transition can be transformed into the first 
order one. However even the picture obtained here 
looks very interesting. 

Let us now discuss some general features of the 
picture obtained focusing on the case of the explicitly 
CP conserving potentials. 

Possible sequences of phase states. 

Thermal evolution can be split into two phases 

I 

(a) EW 

(b) CPc 

(c) CPc EW 



We obviously omitted from this list the sequence 

(EW CPc charged) , which disagrees with to- 
day's Universe. Each possible sequence from the left 
column can be combined with any possible sequence 
from the right column. In the case (b)+(I), the 
history of the Universe contains no phase transition 
at all. 

Regions of parameters allowing for different 
paths of phase evolution. 

We presented a method how to describe regions 
in the space of the zero-temperature parameters of 
the model that lead to each specific type of thermal 
evolution of the Universe (see the shaded areas in 
Figs. EJ. To cast them into the corresponding 



whose properties are rather decoupled from each 
other. 

At the extremely high temperature the vacuum 
state can be either EW symmetric or EW-violating, 
but CP-conserving (CPc). In the former case, which 
we see as the most probable, the Universe during its 
cooling down passes to a CPc phase (EWSB phase 
transition at T — Tew, Fig. left plot). In the sec- 
ond case, cooling down either keeps Universe in the 
CPc phase for a long time or brings it into and then 
out of the EW symmetric phase, that is, first restor- 
ing and then breaking again the EW symmetry. 

As the temperature goes down, the system enters 
the lower-temperature stage of its evolution, which 
in our analysis corresponds to a transition from the 
second to the first sheet of the (/z, 5) plane. The 
mass parameters of the system keep changing, and 
the phase state of system can either vary continuously 
or exhibit one or two additional phase transitions. In 
these transitions our system can either evolve through 
a charged vacuum phase (with very unusual proper- 
ties), or cross the sCPv phase, or stop in the sCPv 
phase state (in this case one should observes spon- 
taneous CP-violation in the Higgs sector in today's 
world) or finally cross the line of the fc-symmetry of 
the CPc phases. In our approximation only the tran- 
sition CPc— >CPc through the fc-symmetric state (at 
k 7^ 1) is of the first order. We therefore list all pos- 
sible sequences of phase states indicating in all cases 
the type of phase transition: 



(I) 
(II) 

(III) (43) 
(IV) 
(V) 

I 

regions of observables, such as masses and couplings 
constants, is a natural task for continuation of this 
work. 

Rearrangement of particle mass spectrum. 

In most examples considered here the value of tan (3 
changes strongly. It jumps at the first order phase 
transition and shows a continuous but very sizable 
evolution for the other cases. In most of our examples 
we saw tan (3 increase towards high temperatures, i.e. 
in the past, but a simple change of the basic Higgs 
fields 1 <-* 2 leads to tan/3 — * 1/ tan/3, that is, to de- 
creasing tan j3. In the latter case equations (TT5)) shows 
that the fermion mass spectrum within one genera- 
tion can be rearranged. Under these circumstances, 



CPc 

CPc — > charged —> CPc 
> ^ CPc-^CPc2 
CPc ^ sCPv 
CPc ±L> sCPv ^ CPc 



18 



it is possible that in the past the decay t — > Wb was 
suppressed, and W — > tb decay was allowed or even 
that the 6-quark was heavier than i-quark. 

In a sequence of several phase transitions, it is only 
the first one, the EWSB, that must take place at the 
electroweak temperature scale. It might happen that 
the other phase transitions take place at much lower 
temperatures. It means that there exists a possi- 
bility that the last phase transition took place rel- 
atively lately in the history of the Universe. If this 
is the case, then the possible rearrangement of the 
quark mass spectrum could have even more spectac- 
ular effects. For example, if the Universe lived long 
enough in an intermediate phase with < m u , 
then the proton could be lighter than the neutron 
and could even decay into it during this intermedi- 
ate stage, which has profound cosmological conse- 
quences. More delicate, even weak variation of rela- 
tion between hadron masses at low temperatures (for 
example, within some stars) can influence for process 
of nucleosynthesis. 

Another interesting opportunity, which can be re- 
alized in many cases, is a non-monotonic dependence 
of masses of particle on temperature, when they start 
from zero at EWSB phase transition, grow and over- 
shoot their today's values and drop down after sub- 
sequent phase transitions. 

In addition, the rearrangement of fermion masses 
can be viewed as yet another phase transition in 
fermion subsystem, with own fluctuations, etc. The 
study of this possible phase transition goes beyond 
the approach developed in this paper. 

Possible relations to cosmology. 

1. Different phenomena discussed here can give rise 
to new effects in the structure of Cosmic Microwave 
Background radiation and other cosmological observ- 
ables. Feasibility of their observation is a subject for 
future studies. 

The cases with new phase transitions in addition 
to the standard EWSB lead to additional stages in 
the early history of the Universe with strong fluctua- 
tions near the phase transition points. For example, 
in many cases, see Figures [5j there exists either 



a meta-stable local minimum state or other extrema 
just above the vacuum state. Possible virtual tran- 
sitions to these states can enhance fluctuations and 
their observable effects. 

2. If the charge- breaking vacuum state was indeed 
an intermediate stage of the evolution of the Uni- 
verse, then a number of unexpected effects appear 
and they can strongly influence the modern situa- 
tion. First of all, in the charge-breaking phase all 
the gauge bosons are massive, and electric charge is 
not conserved. Also, in this phase the long ranged 
forces, like electromagnetism in our world, are absent, 
and the only long range force is gravity. So, the lo- 
cal electric neutrality of medium can be strongly vio- 
lated. After a phase transition to the modern charge- 
conserving vacuum, one could have strong deviations 
from the average electroneutrality, originating from 
the charge-breaking vacuum phase. These are not 
the standard charge fluctuations, and they will result 
in a strong relative motion of separate parts of the 
Universe, which can result either in strong mixing 
and averaging of matter or in production of struc- 
tures like caustics (proto-galaxies) . The restoration 
of the electric neutrality can go on during a very long 
time after the phase transition to neutral vacuum. 

3. In the standard approach the temperature of 
the phase transition is unavoidably set by the elec- 
troweak scale. In our model the same is valid for the 
first EWSB transition. However, the temperature 
of the last phase transition can be sufficiently low, 
that is, the "zero point" in Figs. [5J can be close to 
the phase separation line. Certainly, for a detailed 
description of such situation our approximation must 
be improved. 
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APPENDIX A: TEMPERATURE CORRECTIONS TO POTENTIAL FOR THE CASE OF HARD 

VIOLATION OF Z 2 SYMMETRY 



In the general case of hard violation of Z 2 symmetry all mass terms have non-zero corrections. In addition 
to the usual one loop diagrams of form of Fig. [IIJ one adds diagrams with a specific bilinear vertex generated 
by the x term. We present the resulting mass corrections for the case of explicitly CP-conserving potential 
assuming field renormalization cf>i — > tfii/y/l — x 2 , taking into account a number of x vertex insertions in each 
line. Instead of eq. (TT2)) we obtain after simple algebra 

2d = 3Ai + 2A 3 + A 4 + 6x(2X 6 + A 7 ) + x 2 (6X 3 + 3A 4 + 3A 2 ) + 6k 3 \ 7 , 
2c 2 = 3A 2 + 2A 3 + A 4 + 6x(A 6 + 2A 7 ) + x 2 {6\ 3 + 3A 4 + 3A 2 ) + 6x 3 A 6 , (Al) 
2ci 2 = [6(A 6 + A 7 ) + x{\ 3 + 2A 4 + A 5 )](l + x 2 ) + xr(3Ai + 2A 2 + 4A 3 + 2A 4 ) . 

Here we include the diagram factor in definition of x. 

These equations show similar to that it was shown in sect.Hllthat the curve of physical states is now a straight 
ray in the in the 4-dimensional space {m\ 1 , m 22 , Rem 2 2 , Imm 2 2 ) and consequently in the 3-dimensional space 
(Refi , Im^i, S) (definition of quantities i?e/i, Tro/i is evident in analogy with (|7b[) ). 



APPENDIX B: DECOMPOSITION OF THE POTENTIAL AROUND ARBITRARY EWV 

EXTREMUM 

We present here useful equation Q for decomposition of the potential around arbitrary EWv extremum N, 
obtained with the aid of eq. ([^1) . It is valid for the most general potential (J3J)- The complexity of the potential 
in this form is concentrated in the values As ; 6,7 and y 3 . 

V = £ c N xt + Vi( Xi - y hN ) + K ■ V{4>, N) ; 
V{<t>, N)=x 1 y 2 +x 2 yi -x 3 y& -2/3^3t = 

= (0l(02>Ar-02('/'l)7v) t (0l('/'2)w-02(0l)Af); ^- B1 ^ 

1Z = for charged extremum 

for neutral extremum N . 



K = 2± 



yi + V2 



N 



Here M„± is the squared mass of charged Higgs boson (it can be negative if considered extremum is not 
minimum of potential). The quantity T)(<fr, Af) can be treated as the the "distance" between some set of fields 
and extremum N with (4>i) = {4>)iN ■ By construction (T>) > for any classical values of </>j. 
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APPENDIX C: PARAMETERS FOR FIGURES 

All the figures, showing evolution of physical parameters are obtained by numerical calculations. Here we 
present the parameres, used for the calculations: 



Figure 


Ai 


A 2 


A 3 


A 4 


A 5 
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Fig. [3] 


0.076 


1.22 


1.69 


0.275 


-1.235 


0.39 


-0.17 


Fig. a 


0.222 


3.556 


1.11 


0.205 


0.125 


0.12 


0.02 


Fig. E 


0.488 


2.475 


0.9 


-0.635 


-0.125 


0.175 


-0.17 


Fig. E] 


0.222 


3.556 


1.12 


0.208 


0.428 


0.12 


0.014 


Fig. E| 


0.197 


3.16 


1.21 


0.6 


1.08 


0.065 


-0.05 


Fig.nni 


0.22 


3.556 


1.12 


0.218 


0.418 


0.15 


0.014 


Fig.[[J 


0.217 


3.48 


1.13 


6.68 


-5.32 


0.145 


-0.18 



